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Critical phenomena and quantum phase tran- 
sitions are paradigmatic concepts in modern con- 
densed matter physics. A central example in 
the field of mesoscopic physics is the localization- 
delocalization (metal-insulator) quantum phase 
transition driven by disorder — the Anderson 
transition. -"^ Although the notion of localization 
has appeared half a century ago, this field is 
still full of surprising new developments. The 
most recent arenas where novel peculiar localiza- 
tion phenomena have been studied are graphene^ 
and topological insulators, ^^s^i^i^ii i.e., bulk insu- 
lators with delocalized (topologically protected) 
states on their surface. Besides exciting phys- 
ical properties, the topological protection ren- 
ders such systems promising candidates for a va- 
riety of prospective electronic and spintronic de- 
vices. It is thus of crucial importance to un- 
derstand properties of boundary metallic modes 
in the realistic systems when both disorder and 
interaction are present. Here we find a novel 
critical state which emerges in the bulk of two- 
dimensional quantum spin Hall (QSH) systems 
and on the surface of three-dimensional topolog- 
ical insulators with strong spin-orbit interaction 
due to the interplay of nontrivial Z2 topology and 
the Coulomb repulsion. At low temperatures, 
this state possesses a universal value of electri- 
cal conductivity. In particular, we predict that 
the direct QSH phase transition occurs via this 
novel state. Remarkably, the interaction-induced 
critical state emerges on the surface of a three- 
dimensional topological insulator without any ad- 
justable parameters. This "self-organized quan- 
tum criticality" is a novel concept in the field of 
interacting disordered systems. 

The universality of critical phenomena has been stud- 
ied both in the area of strongly correlated systems and 
in mesoscopics. It is now established that disordered 
electronic systems can be classified into 10 symmetry 
classes (for review see Ref. [l|) . Very generally, the local- 
ization properties are determined by the symmetry class 
and dimensionality of the system. The critical behav- 
ior of a system depends also on the underlying topol- 
ogy. This is particularly relevant for topological insula- 
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The famous example of a topological insulator is a 



two-dimensional (2D) system on one of quantum Hall 
plateaus in the integer quantum Hall effect (QHE). Such 
a system is characterized by an integer (Chern number) 
n = ...,—2,-1,0,1,2,... which counts the edge states. 
Here the sign determines the direction of chiral edge 
modes. The integer quantum Hall edge is thus a topo- 
logically protected ID conductor realizing the group Z. 

Another (Z2) class of topological insulators24,iS, can 
be realized in systems with strong spin-orbit interac- 
tion and without magnetic field — and was discovered 
in HgTe/HgCdTe structures in Ref. (see also Ref. 8). 
Such systems were found to possess two distinct insu- 
lating phases, both having a gap in the bulk electron 
spectrum but differing by the edge properties. While the 
normal insulating phase has no edge states, the topolog- 
ically nontrivial insulator is characterized by a pair of 
mutually time-reversed delocalized edge states penetrat- 
ing the bulk gap. Such state shows the quantum spin 
Hall (QSH) effect which was theoretically predicted in a 
model system of graphene with spin-orbit coupling4iii^ 
The transition between the two topologically nonequiv- 
alent phases (ordinary and QSH insulators) is driven by 
inverting the band gap?^ The Z2 topological order is ro- 
bust with respect to disorder: since the time-reversal in- 
variance forbids backscattering of the edge states at the 
boundary of QSH insulators, these states are topologi- 
cally protected from localization. 

A related three-dimensional (3D) Z2 topological insu- 
lator was discovered in Ref. where crystals of Bii_a;Sba; 
were investigated. The boundary in this case gives rise 
to a 2D topologically protected metal. Similarly to 
2D topological insulators, the inversion of the 3D band 
gap induces an odd number of the surface 2D modesJ^ 
These states in BiSb have been studied experimentally 
in Refs. 't"^ and [T^. Other examples of 3D topological 
insulators include BiTe and BiSe systemsJii^iii^ 

Both in 2D and 3D, Z2 topological insulators are band 
insulators with the following properties: (i) time rever- 
sal invariance is preserved (unlike ordinary quantum Hall 
systems); (ii) there exists a topological invariant, which 
is similar to the Chern number in QHE; (iii) this invari- 
ant belongs to the group Z2 and reflects the presence or 
absence of delocalized edge modes (Kramers pairs) 4 For 
the sake of completeness we overview the full classifica- 
tion of topological insulators and superconductorsiS in 
Supplementary Information. 

In this paper, we consider the effect of interactions 
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FIG. 1: Schematic scaling functions for the conductivity of 2D disordered systems of symplectic symmetry class. The plotted 
beta functions l3{g) — dg/dlnL determine the flow of the dimensionless conductivity g with increasing system size L (as 
indicated by the arrows). The upper two panels show the beta functions for ordinary spin-orbit systems which are not 
topologically protected (left: no interaction; right: Coulomb interaction included). The lower two panels demonstrate the 
scaling for topologically protected Dirac fermions (left: no interaction; right: Coulomb interaction included). 



on Z2 topological insulators belonging to the symplec- 
tic symmetry class, characteristic to systems with strong 
spin-orbit interaction. In particular, we predict a novel 
critical state which emerges in a 2D system due to the 
interplay of nontrivial topology and the Coulomb inter- 
action. 

Let us start with reviewing the localization properties 
of 2D systems of symplectic symmetry class All without 
Coulomb interaction. In conventional spin-orbit systems 
(e.g. semiconductors with spin-orbit scattering), there 
are two phases: metal and insulator with the Ander- 
son transition between them, see Fig. [TJ A qualitatively 
different situation occurs in a single species of massless 
Dirac fermions in a random scalar potential. This sys- 
tem also belongs to the symplectic symmetry class but its 
metallic phase is "topologically protected" whatever dis- 
order strength. This state has been recently predicted 
for disordered graphene with no spin- and no valley- 
mixingiiiii^ii^ The absence of localization in this model 
as well as the "supermetal" scaling (Fig. [1]) have been 
confirmed in numerical simulations Although a gen- 
uine single Dirac fermion cannot be realized in a truly 
2D microscopic theory because of the famous "fermion 
doubling" problem, it emerges on the surface of a 3D 
topological insulator. 

The 3D topological insulators are characterized by the 
inverted sign of the gap (band inversion). This generates 
the surface states, as was first pointed out in Ref. IH. 
The effective 2D surface Hamiltonian has a Rashba form 
(see the derivation in Supplementary Information) and 
describes a single species of 2D massless Dirac particles 



(cf. Ref. [23I) . It is thus analogous to the Hamiltonian of 
graphene with just a single valley. In the absence of in- 
teraction, the conductivity of the disordered surface of a 
3D topological insulator therefore scales to infinity with 
increasing the system size. This behavior defines the 
topologically-protected "supermetallic phase" discussed 
above. 

Let us now "turn on" the Coulomb interaction between 
electrons. Since a topological insulator is characterized 
by the presence of non-localizable surface states, its ro- 
bustness with respect to interaction can be reformulated 
in terms of the absence of localization of the boundary 
states. At this point it is worth recalling the celebrated 
example of a 2D topological insulator, the QHE insula- 
tor, in which the interaction cannot destroy the chiral ID 
edge modes on the boundary of a 2D QHE sample. Fur- 
thermore, the two consequent QHE topological insulators 
(QHE plateaus) are separated by a delocalized (critical) 
2D state. Since the topological insulator phases are ro- 
bust, the interaction is not capable of localizing electrons 
in this state separating the two topological insulators. 

Having in mind this analogy with the QHE, one can 
expect that the delocalized states in Z2 topological insu- 
lators are not destroyed by the interaction either. Indeed, 
arguments in favor of the stability of Z2 topological insu- 
lators with respect to interactions were given in Refs. |3 
and [23. Below we demonstrate the Z2 topological or- 
der and discuss its implications in 2D and 3D interacting 
systems in the presence of disorder. 

We first consider the interacting massless Dirac elec- 
trons on the surface of a 3D topological insulator. With- 
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FIG. 2: Left: schematic illustration of the hollow cylinder sample used for proving the absence of localization on the surface of 
a 3D topological insulator (see Supplementary Information). Right: the energy spectra of a clean 2D system on the surface of 
the cylinder with zero flux and with half of the magnetic flux quantum penetrating the cylinder. 



out interaction, the surface states are delocalized in the 
presence of arbitrarily strong potential disorder. In Sup- 
plementary Information we demonstrate that the interac- 
tion cannot fully localize the surface Dirac fermions. This 
is achieved by considering the topological insulator of a 
hollow cylinder geometry threaded by half of the mag- 
netic flux quantum, see Fig. O Our 2D problem then 
reduces to the quasi-lD model with an odd number of 
channels. Full 2D localization would be in contradiction 
with known results on the absence of localization in such 
quasi-lD symplectic wires. Since delocalization in quasi- 
lD geometry survives in the presence of interaction, this 
is also true for the 2D interacting Dirac electrons on the 
surface of a 3D topological insulator. 

Can the topologically protected 2D state be a super- 
metal {g oo) as in the noninteracting case? To answer 
this question we employ the perturbative renormalization 
group (RG) applicable for large conductivity g 1. 

It is well known that in a 2D diffusive system the inter- 
action leads to logarithmic corrections to the conductiv- 
ity.— These corrections (together with the interference- 
induced ones) can be summed up with the use of RG tech- 
niqucj^i^liS^ The one-loop equation for renormalization 
of the conductivity in the symplectic class with Coulomb 
interaction has the following form: 

dg 
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d\nL 

where N is the number of degenerate species ( "flavors" : 
spin, isospin, . . . ) and L is the system size. 

The first term, N/2, describes the effect of weak an- 
tilocalization due to disorder (this term exists also in the 
absence of interaction) for N parallel conductors. The 
second term, —1, is induced by the Coulomb interaction 
in the singlet channel and has a localizing effect: it sup- 
presses the conductivity. The last term on the r.h.s. of 
Eq. ([T]) is due to the interaction in the multiplet (in the 
flavor space) channel. This term yields a positive (antilo- 
calizing) correction to the conductivity. The multiplet 
interaction parameter J-' is itself subject to renormaliza- 



In the degenerate case N > 1 (as, for example, in 
graphene with no spin- and valley- mixing where TV = 4), 
the beta function ([T]) is positive corresponding to the "su- 
permetal" phase. The situation is essentially different for 
2D states on the surface of a 3D Z2 topological insulator 
where we have a symplectic system with = 1. Accord- 
ing to Eq. ID), the negative interaction- induced term in 
f3{g) now dominates while the multiplet term is absent. 
Therefore, for g ^ 1 the conductance decreases upon 
renormalization. This means that due to interaction the 
supermetal fixed point becomes repulsive. 

Thus, on one hand, we encounter the tendency to lo- 
calization due to the interaction. On the other hand, 
the states on the surface of the topological insulator are 
topologically protected from the localization. At g ^ 1 
the topological protection reverses the sign of the /3 func- 
tion, similarly to the ordinary QHE. As a result, a critical 
point emerges due to the combined effect of interaction 
and topology, see Fig.[T] This type of criticality should be 
contrasted with the QHE criticality which exists already 
without interactions. In our case, even in the absence 
of the critical state in a non-interacting model, the criti- 
cality is inevitably established in the realistic interacting 
systems. This novel interaction-induced critical state is 
the major result of our paper. 

Let us now return to the case of 2D Z2 topological 
insulators. As discussed above, without interaction dis- 
order was found to induce a metallic phase separating 
the two (QSH and ordinary) insulating phases. The 
conductance in the metallic phase scales to infinity be- 
cause of the weak antilocalization specific to the symplec- 
tic symmetry class. The schematic phase diagram for the 
noninteracting case is shown in Fig. [3] (left panel). 

What would change in this phase diagram when inter- 
action is taken into account? The answer follows from 
Eq. The 2D disordered QSH system contains only a 
single flavor of particles, N = 1. Indeed, the spin-orbit 
interaction breaks the spin-rotational symmetry, whereas 
the valleys are mixed by disorder. As a result, the super- 
metal phase does not survive in the presence of Coulomb 
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band gap band gap 

FIG. 3: The phase diagrams of a disordered 2D system demonstrating the QSH effect. Left: noninteracting case. In contrast to 
the clean case, the two topologically nonequivalent phases (QSH and ordinary insulators) have been found to be separated by 
the metallic phase.— The transition'^ between metal and any of the two insulators occurs at the critical value of dimensionless 
(in units /h) conductivity g = g* ~ 1.4; both transitions are believed to belong to the same universality class. For g < g* 
all bulk states are eventually localized in the limit of large system, while for g > g* antilocalization drives the system to the 
"supermetallic" state, g —> oo. Right: interacting case (with Coulomb interaction not screened by external gates). Interaction 
"kills" the supermetallic phase. As a result, the two insulating phases are separated by the critical line. 



interaction: at g 3> 1 the interaction-induced localization 
wins. This is analogous to the case of the surface of a 3D 
topological insulators discussed above. 

The edge of a 2D topological insulator is protected 
from the full localization, as was discussed already in the 
pioneering works by Kane and Melej^,)! In the presence of 
interaction, the counter-propagating edge modes consti- 
tute the Luttinger liquid. As disorder-induced backscat- 
tering in this liquid is forbidden by the time-reversal 
symmetry, the dimensionless conductance of this ID sys- 
tem is equal to unity. This means that the topologi- 
cal distinction between the two insulating phases (ordi- 
nary and QSH insulator) is not destroyed by the inter- 
action, whereas the supermetallic phase separating them 
disappears. Therefore we conclude that the transition 
between two insulators occurs through an interaction- 
induced critical state, see Fig. [3] (right panel). 

We have thus found that in the two types of systems 
with strong spin-orbit coupling the Coulomb interaction 
induces novel 2D critical states. This happens, first, at 
the boundary of 3D topological insulators and, second, in 
the bulk of 2D QSH systems, where the critical state sep- 
arates the two topologically distinct insulating phases. In 
the first case the system can be described by a 2D inter- 
acting symplectic sigma-model with the Z2 topological 
term. The two critical states have much in common: (i) 
symplectic symmetry, (ii) Z2 topological protection, (iii) 
interaction-induced criticality, and (iv) conductivity of 
order unity (probably universal). This suggests that the 
corresponding fixed points might be equivalent. The ex- 
istence of the proposed critical states can be verified in 
transport experiments on the semiconductor structures 
with possible gap inversion. 
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valuable discussions. The work was supported by 
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the EUROHORCS/ESF (IVG), and by Rosnauka grant 
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I. SUPPLEMENTARY INFORMATION 

Classification of topological insulators 

The full classification (periodic table) of topological 
insulators and superconductors for all te n sy mmetry 
classes-^'-- ^g^g developed in Refs. [H and [H. In this 
Section we overview the classification of topological in- 
sulators closely following Refs. 'si' and'ss' and discuss the 
connection between the classification schemes of these pa- 
pers. Very generally, the classification of topological in- 
sulators in dimension D can be constructed by studying 
the Anderson localization problem in a D-1 disordered 
system: absence of localization of surface states implies 
the topological character of the insulator)^ 

All symmetry classes of disordered systems^i^^ can be 
divided into two groups: {A,AIII} and {all other}. The 
classes of the big group are labeled by p = 0, 1, . . . , 7. 
Each class is characterized by (i) Hamiltonian symmetry 
class Hp; (ii) symmetry class of the classifying space used 
by Kitaev"^^ Rp; (iii) symmetry class Sp of the compact 
sector A4.P oi the sigma-model manifold. The symmetry 
class Rp of the classifying space of reduced Hamiltoni- 
ans characterizes the space of matrices obtained from the 
Hamiltonian by keeping all eigenvectors and replacing all 
positive eigenvalues by -1-1 and all negative by —1. Note 
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TABLE I: Symmetry classes and "Periodic Table" of topological insulators; - The first column enumerates the symmetry 
classes of disordered systems which are defined as the symmetry classes Hp of the Ifamiltonians (second column) . The third 
column lists the symmetry classes of the classifying spaces (spaces of reduced Hamiltonians) .— The fourth column represents 
the symmetry classes of a compact sector of the sigma-model manifold. The fifth column displays the zeroth homotopy group 
T^oiRp) of the classifying space. The last four columns show the possibility of existence of Z and Z2 topological insulators in 
each symmetry class in dimensions d = 1, 2, 3, 4. 
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Here and below cyclic definition of indices {0,1,. 
(mod 8) and {0', 1'} (mod 2) is assumed. 

For the classification of topological insulators it is im- 
portant to know homotopy groups tt^ for all symmetry 
classes. In Table U we list TTo{Rp); other tt^ are given by 



T^diRp) — T^oiRp+d) 



(3) 



The homotopy groups Tr^ have periodicity 8 (Bott peri- 
odicity). 

There are two ways to detect topological insulators: by 
inspecting the topology of (i) classifying space Rp or of 
(ii) the sigma-model space Sp. 

(i) Existence of topological insulator (TI) of class p in 
d dimensions is established by the homotopy group ttq 
for the classifying space Rp-d- 



TI of the type Z 
TI of the type Z2 



T^oiRp-d) 



Z2 



(4) 



(ii) Alternatively, the existence of topological insulator 
of symmetry class p in d dimensions can be inferred from 
the homotopy groups of the sigma-model manifolds, as 
follows: 



fTI of the type Z 
I TI of the type Z2 



T^d{Sp) = Z 

''^d-iiSp) = Z2 



(5) 



This criterium is obtained if one requires existence of 
"non-localizable" boundary excitations. This may be 
guaranteed by either Wess-Zumino term in d — 1 dimen- 
sions [which is equivalent to the Z topological term in d 



dimensions, i.e. Trd{Sp) = Z] for a QHE-type topological 
insulator, or by the Z2 topological term in d — 1 dimen- 
sions [i.e. TTd-i{Sp) = Z2] for a QSH-type topological 
insulator. 

The above criteria (i) and (ii) are equivalent, since 



and 



TTdiSp) = -KdiRi-p) = 7ro(i?4-p+d)- 



7ro(i?p) 




(6) 



(7) 



In this paper we focus on 2D systems of symplectic 
(All) symmetry class. One sees that this is the only 
symmetry class out of ten classes that supports the ex- 
istence of Z2 topological insulators both in 2D and 3D. 
The effect of interaction on Z2 topological insulators and 
superconductors in classes DIII (in 2D) and CII (in 3D) 
will be considered elsewhere. 

Surface states of 3D topological insulators 

In this section we consider the surface of 3D topological 
insulators and the derive the effective surface Hamilto- 
nian. The realistic microscopic Hamiltonian of such sys- 
tems with strong spin-orbit interaction can be modelled 
by the 3D massive Dirac Hamiltonian, see, e.g., Refs. [isl 
and [1^. We start with the general form of a 3D Dirac 
Hamiltonian 



H3D — 



-M o-p' 
crp M . 



(8) 



This Hamiltonian is a matrix in the 4x4 space formed 
by the spin and pseudospin (sublattice) spaces. The in- 
terface between the semiconductor and the vacuum is 
described by sending the mass M to infinity. 
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Let us consider a flat interface at x = 0. The edge 
state with zero energy decays into the bulk {x < 0): 

VP = e''- (^^ . (9) 

Acting by the Hamiltonian ([S]) on this wavefunction, we 
obtain the relation between the two components of the 
spinor, x — io'x'>p- For an arbitrary surface characterized 
by the normal vector n, the general boundary condition 
reads 

X = i<jnip. (10) 

From Eqs. ([5]) and ([TO)) we obtain the effective 2D surface 
Hamiltonian in the form of the Rashba Hamiltonian in a 
curved space 

^/"surf = ^ [ctp, crn] = ^ + i (n[p X cr] + [p X cr]n) . 

(11) 

This 2D Hamiltonian describes a single species of 2D 
massless Dirac particles and is thus analogous to the 
Hamiltonian of graphene with just a single valley. 

Absence of localizations of 2D surface states 

Consider the 2D system formed at the surface of a 3D 
topological insulator and described by the Hamiltonian 
PTjl . Let us first prove the absence of localization in the 
non-interacting case. Assume that all the states in 2D 
are localized with some localization length ^. Consider a 
hollow cylinder with all dimensions much larger than ^ 
pierced by the Aharonov-Bohm magnetic flux $ = hc/2e 
(half of the flux quantum). This value of $ does not 
break the time-reversal symmetry leaving the system in 
the symplectic class. In the absence of disorder we can 
characterize the surface states by the momentum k along 
the cylinder axis and by the integer angular momentum 
n. The energy of such a state is given by 

E = y^k^ + in/R)'^. (12) 



The channels with positive and negative n are degen- 
erate, while n — channel is not. Thus the cylinder 
sustains an odd number of conducting channels both on 
the inner and outer surface at any value of chemical po- 
tential. 

Let us now include disorder and show the absence of lo- 
calization in quasi-one-dimensional (qlD) symplectic sys- 
tem with odd number of channels. The scattering matrix 
of such a qlD wire has the form 

with transmission and reflection amplitudes as its entries. 
The blocks r and r' are square matrices of the size deter- 
mined by the number of channels. Time-reversal symme- 
try of the symplectic type imposes the following restric- 
tions on the amplitudes entering the matrix S: 

r = -r'^, r' = -r'^, t ^ t'^ . (14) 
Calculating the determinant of the both sides of the flrst 
identity and taking into account the odd size of the ma- 
trix r, we obtain det r = 0. This implies a zero eigenvalue 
of r and hence the existence of a channel with perfect 
transmission. We conclude: in a qlD wire of symplectic 
symmetry with an odd number of channels one channel 
always remains delocalized.— i^i^i^ 

Applying the qlD result to the cyhnder constructed 
above we immediately come to the controversy: in spite 
of assumed 2D localization on the surface, the infinitely 
long cylinder possesses two (inner and outer) conducting 
channels. This proves the absence of localization in 2D. 

The proof can be generalized to include the Coulomb 
interaction. We assume the temperature to be much 
smaller than the inverse time of electron propagation 
through the system. At such low temperatures the in- 
elastic scattering of electrons is negligible and we can de- 
scribe the transport by the scattering matrix calculated 
at the Fermi energy and accounting for virtual processes. 
The symmetry properties of this S matrix are unchanged 
and hence the above proof applies. 



^ F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 

(2008) . 

^ A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009) . 

^ C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

^ C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

^ B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 

314, 1757 (2006). 
® M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buhmann, 

L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 

766 (2007). 

^ D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, 



and M. Z. Hasan, Nature 452, 970 (2008). 
* A. Roth, C. Briine, H. Buhmann, L. W. Molenkamp, J. 
Maciejko, X.-L. Qi, and S.-C. Zhang, Science 325, 294 
(2009). 

8 D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal , J. H. Dil, 
J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. 
S. Hor, R. J. Cava, and M. Z. Hasan, Science 323, 919 
(2008). 

° P. Roushan, J. Seo, C. V. Parker, Y. S. Hor, D. Hsieh, 

D. Qian, A. Richardella, M. Z. Hasan, R. J. Cava, and A. 

Yazdani, Nature 460, 1106 (2009). 
^ D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, J. 

Osterwalder, L. Patthey, J. G. Checkelsky, N. P. Ong, A. 

V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. 

Cava, and M. Z. Hasan, Nature 460, 1101 (2009). 



7 



12 H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. 

Zhang, Nature Phys. 5, 438 (2009) 

Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S. K. Mo, 

X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. 

Zhang, L R. Fisher, Z. Hussain, and Z.-X. Shen, Science 

325, 178 (2009). 
1" L. Sheng, D. N. Sheng, C. S. Ting, and F. D. M. Haldane, 

Phys. Rev. Lett. 95, 136602 (2005). 
1^ L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
1® A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- 

wig, Phys. Rev. B 78, 195125 (2008). 
1'' P. M. Ostrovsky, L V. Gornyi, and A. D. Mirlin, Phys. 

Rev. Lett. 98, 256801 (2007). 
1* P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin, Eur. Phys. 

J. Special Topics 148, 63 (2007). 
1^ S. Ryu, C. Mudry, H. Obuse, and A. Furusaki, Phys. Rev. 

Lett. 99, 116601 (2007). 
2° J. H. Bardarson, ,J. Tworzydlo, P. W. Brouwcr, and C. W. 

J. Beenakker, Phys. Rev. Lett. 99, 106801 (2007). 

21 K. Nomura, M. Koshino, and S. Ryu, Phys. Rev. Lett. 99, 
146806 (2007). 

22 M. L Dyakonov and A. V. Khaetskii, JETP Lett. 33, 115 
(1981). 

23 B. A. Volkov and O. A. Pankratov, JETP Lett. 42, 178 
(1985). 

2* S.-S. Lee and S. Ryu, Phys. Rev. Lett. 100, 186807 (2008). 



2® B.L. Altshulor and A.G. Aronov, in Electron- electron in- 
teractions in disordered conductors, edited by A.L. Efros 
and M. PoUak (Elsevier, 1985), p. 1. 

2^ A. M. Finkelstein, Sov. Sci. Rev. A. Phys. 14, 1 (1990). 

2'' A. Punnoose and A. M. Finkel'stein, Phys. Rev. Lett. 88, 
016802 (2001). 

2* D. Belitz and T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 
(1994). 

2^ M. Onoda, Y. Avishai, and N. Nagaosa, Phys. Rev. Lett. 

98, 076802 (2007). 
3° H. Obuse, A. Furusaki, S. Ryu, and C. Mudry, Phys. Rev. 

B 76, 075301 (2007). 

31 M. R. Zirnbauer, J. Math. Phys. 37, 4986 (1996). 

32 A. Ahland and M. R. Zirnbauer, Phys. Rev. B 55, 1142 
(1997). 

33 A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- 
wig, AIP Conf. Proc. 1134, 10 (2009). 

3* A. Yu. Kitaev, AIP Conf. Proc. 1134, 22 (2009). 

35 M.R. Zirnbauer, Phys. Rev. Lett. 69, 1584-1587 (1992). 

3^ A.D. Mirlin, A. Miiller-Groeling, and M.R. Zirnbauer, 

Ann. Phys. 236, 325 (1994). 
3^" T. Ando and H. Suzuura, J. Phys. Soc. Jpn. 71, 2753 

(2002). 

3* Y. Takane, J. Phys. Soc. Jpn. 73, 1430 (2004). 



